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Abstract 

An operator that governs the discrete time evolution of the velocity distribution of an out-of-equilibrium ideal gas will be 
presented. This nonlinear map, which conserves the momentum and the energy of the ideal gas, has the Maxwellian Velocity 
Distribution (MVD) as an asymptotic equilibrium. Moreover, the system displays the increasing of the entropy during the decay 
to the MVD. 
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I. Introduction 

^jPj A scheme inspired in economic systems that has recently been proposed [1] to explain the attractivity, and then the ubiquity, 
^ of the exponential (Boltzmann-Gibbs) Distribution (BGD) reads as follows: Let p{m)dm denote the PDF (probability density 
function) of money in a multi-agent economic system, i.e. the probability of finding an agent of the ensemble with money 
between m and m + dm. Consider now the discrete time evolution of an initial money distribution po{m) at each time step 
, n under the action of an operator T, which represents the average effect on the system of many random binary interactions 
^ (with number of the order of the system size) between pairs of agents exchanging their money. Thus, the system evolves from 
.time n to time 71 + 1 to asymptotically reach the equilibrium wealth distribution pf{m), i.e. 

<^ ■ lim r" {po{m)) Pf{x) = BGD = Se-^"" with S =< po . 



,In this case, Pf{m) is the exponential distribution (BGD) with the same average wealth < pf >= S^^ than the initial one 
< Po due to the local and total money conservation [2]. The mathematical properties of operator T have been disclosed in 
,Ref. [3]. Hence, this framework not only puts in evidence that the BGD is the equilibrium distribution if not that in this case 
■ the BGD is asymptotically reached independently of the initial wealth distribution given to the system, a point of view that to 
^ , date was possibly lacking in the literature. 

In this work, we extend this perspective to another problem of the same statistical nature. Our goal is to explain the ubiquity 
of the MVD in ideal gases [4], [5]. In the next section, we explain how to obtain an operator T in the space of velocity 
On distributions in order to explain the decay of any initial velocity distribution to the MVD. Then, the dynamical properties of 
^""^ [ this operator T will be sketched and some examples showing its dissipative behavior are depicted. Finally, our conclusions are 
00 'given. 
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^ \ II. The Model. Nonlinear Map T 

' Consider an ideal gas with particles of unity mass in the three-dimensional (3Z)) space. As long as there is not a privileged 
J> .direction in the equilibrium, we can take any direction in the space and to study the discrete time evolution of the velocity 
distribution in that direction. Let us call this axis U . We can complete a Cartesian system with two additional orthogonal 
axis V,W. If Pn{u)du represents the probability of finding a particle of the gas with velocity component in the direction 
• U comprised between u and u + du at time n, then the probability to have at this time n a particle with a 3-D velocity 
{u,v,w) will be pn{u)pn{v)pn{w) . The particles of the gas collide between them, and after a number of interactions of the 
order of system size, a new velocity distribution is attained at time n + 1. Concerning the interaction of particles with the 
bulk of the gas, we make two simplistic and realistic assumptions in order to obtain the probability of having a velocity x 
in the direction U at time n + 1: (1) Only those particles with an energy bigger than at time n can contribute to this 
velocity x in the direction U, that is, all those particles whose velocities {u,v,w) verify + + > x^; (2) The new 
velocities after collisions are equally distributed in their permitted ranges, that is, particles with velocity (u, v, w) can generate 
maximal velocities ±11 max = iVu^ + v'^ + w^, then the allowed range of velocities [—Umax, Umax] measures 2|[/„iaa;|, and 
the contributing probability of these particles to the velocity x will be Pn{u)pn{v)pn{w) / {2\Umax\)- Taking all together we 
finally get the expression for the evolution operator T. This is: 

P„.t(x) = T,„(x) ^( ( ( ^»(;K(^K(^) d,d.d.. . 
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Let us remark that we have not made any supposition about the type of interactions or colHsions between the particles and, 
in some way, the equivalent of the Boltzmann hypothesis of molecular chaos would be the two simplistic assumptions we 
have stated on the interaction of particles with the bulk of the gas. But now a more clear and understandable framework than 
those usually presented in the literature appears on the scene. In fact, the operator T conserves in time the energy and the null 
momentum of the gas. Moreover, for any initial velocity distribution, the system tends towards its equiUbrium, i.e. towards the 
MVD. This means that 

lim T" {po{x)) -> pf{x) = MVD {ID case) . 

In the presentation, we will pass review to the properties that (until today) have been collected here [6]. Other additional 
properties concerning the increasing of the entropy with time will also be given. Let us recall at this point the result that 
resumes the behavior of the system: 

Conjecture: For any p ^ B, with B the set of functions whose norm | \p\ \ is equal to unity, with finite mean energy, < ,p >, 
and verifying lim„^oo \\T^p{x)~ IJ-i^)]] = 0, the limit ij,{x) is the fixed point = y^e~"^ , with a = (2 < x^,p >)~^. 

In physical terms, it means that for any initial velocity distribution of the gas, it decays to the Maxwellian distribution, which 
is just the fixed point of the dynamics. Recalling that < x^,p >= fcr, with k the Boltzmann constant and r the temperature 
of the gas, and introducing the mass m of the particles, let us observe that the MVD is recovered in its 3-D format: 
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MVD = Pa(uK(wKH - (^) ' +^'+'"') w^th a = {2kT)-\ 

Some examples showing this behavior are depicted in the Figures 1-2. 
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Fig. 1. (a) p{x) and fJ.{x), (b) Tp{x) and 




(a) (b) 



Fig. 2. (a) p{x) and ^i(x), (b) Tp{x) and fJ,{x). 

III. Conclusions 

In this work, a nonlinear map acting on the velocity distribution space of ideal gases, which gives account of the decay of 
an out-of-equilibrium velocity distribution toward the Maxwellian distribution, has been presented. Some properties and figures 
concerning its dynamical behavior have also been shown. 
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